Spin-current in generic hybrid structures due to interfacial spin-orbit scattering 
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We demonstrate a general principle that hybrid structures of any sort inevitably will give rise to a pure spin- 
current flowing parallel to the interface region when a charge-current is injected. This stems from the broken 
mirror symmetry near the interface which gives rise to spin-orbit coupling that deflects incoming electrons in 
a spin-discriminating fashion. We establish a general analytical condition for the appearance of this effect, 
and calculate the transverse spin-current explicitly using two different models. In addition, we investigate how 
the process of Andreev-reflection influences this phenomenon in the scenario where one of the materials is 
superconducting. 
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Whereas the charge of electrons conventionally has been 
utilized in the field of electronics, the transport of electron 
spin and its belonging detection and manipulation has devel- 
oped into a major research field in modern condensed matter 
physics over the last two decades [ 1 1. Besides its allure from 
a fundamental physics point of view, practical applications re- 
lated to spintronics are already well-established e.g. in the 
commercial harddisk-drive industry. More recently, state-of- 
the-art spintronics revolves around manipulation of domain 
walls in nanowires which could lead to the development of 
so-called race track memory devices Q. 

A central concept in the field of spintronics is the spin- 
current: a flow of spin angular momentum carried by elec- 
trons. It is desirable to find experimentally feasible methods 
to generate spin-currents which are tunable and possible to 
exert control over. In fact, the spin-current is a key part of un- 
derstanding phenomena such as the giant magnetoresistance 
effect |H and current-induced spin-transfer torque |5] |6] . 
One of the main challenges in spintronics is to create pure 
spin-currents for the purpose of spin injection. There exists 
several proposals which revolve around the spin Hall effect, 
spin-current generation transverse to the applied electric field 
due to spin-orbit coupling, to achieve this goal iTT HTTjll . 

The bulk properties of a material are of high importance 
with regard to how spin-polarized transport may be obtained, 
often necessitating the use of materials with intrinsic fer- 
romagnetic correlations in order to spin-polarize a current. 
However, the interface properties in hybrid structures are also 
an essential ingredient in this endevour. In fact, a generic 
property of any hybrid structure is that inversion symmetry is 
broken near the interface: the system is no longer centrosym- 
metric. As a direct consequence of this, a gradient of the elec- 
tric potential will inevitably give rise to antisymmetric spin- 
orbit coupling effects for electrons moving in this region. This 
observation suggests a striking effect. In this Letter, we will 
demonstrate that the broken centrosymmetricity inherit to all 
hybrid structures invariably will generate a pure spin-current 
flowing parallel to the interface even in the absence of any 
magnetic elements or correlations, upon injection of charge- 
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FIG. 1: (Color online) Basic mechanism for the generic spin-current 
in hybrid structures. Centrosymmetricity is broken near the interface 
region, giving rise to an effective electric field. Incoming electrons 
thus feel an effective ^-dependent magnetic field which couples to 
their spin. In this way, a pure transverse spin-current is generated 
near the interface region. In the calculations, we restrict our attention 
to a planar (xy) structure to obtain more transparent analytical results, 
although they may be generalized to 3D structures. 



current {e.g. by current- or voltage-bias). 

Let us explain the physical mechanism which gives rise to 
the general appearance of a transverse pure spin-current in a 
hybrid structure. Consider the model shown in Fig. [TJ a bi- 
layer setup with two different materials separated by an in- 
terface region. Near the interface, mirror symmetry is bro- 
ken and the system is non-centrosymmetric. As a result, an 
incoming electron with velocity v experiences a net electric 
field E due to the asymmetric crystal-field potential, induced 
perpendicularly to the interface. In the rest-frame of the elec- 
tron, an effective magnetic field B^vxE appears. In turn, 
this gives rise to a spin-orbit coupling interaction of the form 
B ■ S °< h ■ (k x a), where h is an interface normal vector. As a 
consequence, the total interface potential matrix V = Vo + Vsoc 
in spin space may be written as a sum of non-magnetic part 
(Vo) and spin-orbit coupling contribution (Vsoc)- This interface 
potential breaks spin symmetry while a transverse charge cur- 
rent is absent due to time-reversal symmetry. Therefore, the 
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injected current generates a pure transverse spin-current as 
a result of the broken centrosymmetricity, which is a generic 
feature of any hybrid structure. 

To illustrate this explicitly, we consider the model shown in 
Fig. [TJsuch that h || x and demonstrate that a pure spin-current 
will flow parallel to the interface, whereas there is no net spin 
transport across the junction. We will show that this transverse 
spin-current arises solely due to the broken centrosymmetric- 
ity, without any magnetic elements incorporated. To begin 
with, we analyze a prototype system consisting of a bilayer 
junction with two normal metals. We assume that there is a 
small region near the interface where spin-orbit coupling ef- 
fects are present due to the lack of inversion symmetry. This 
barrier region is taken to have a finite width L and a potential 
Uq in addition to spin-orbit coupling U soc . In this case, we may 
write: V = [UqI + U soc h ■ (a x k)/k F ]&(x)&(x - L), where 
®(x) is the Heaviside step-function. The wavefunction in the 
left normal metal region, where we will evaluate the gener- 
ated spin-current, reads: = (8 O |,8 a |)e 1 ' ?;c + r°(l,0)e _1< ' A ' + 

r°(Q, \)eT lqx . Here, q = kfeosQ is the wavevector in the nor- 
mal region and kp is the magnitude of the Fermi momen- 
tum. The boundary conditions then demand = \\fM and 
v^l = vJVm at x = 0, whereas y R = y M and y^\\t R = vJVm 
at x = L. Here, and \\i R are wavefunctions in the mid- 
dle and right regions, respectively. The velocity operators 
are given by v_ x = dH/dk x . The above set of equations may 
be solved analytically to provide the scattering coefficients of 
the system, which in turn can be used to calculate the spin- 
currents in the system. The spin-current is in general a tensor, 
with a direction of flow in real space and a polarization in spin 
space. For concreteness, we consider the z-component of the 
polarization and focus on its flow in the x- and y-directions 
which correspond to perpendicular and parallel to the barrier, 
respectively. The general definition reads: 




Inserting the wavefunction \\tp into Eq. ([TJ, we 
note that J x can be written generally as 11241 : J x — 
Jo /2 ddk F cos0/(2m)£ ao [|r°(a0)| 2 - \rf(aB)\ 2 ]. Therefore, 
the general condition for the absence of a spin-current flow- 
ing across the barrier reads: 

£[|^(a0)| 2 - |rj(a9)| 2 ] = with a = ±1, a=t4 • (2) 

aa 

Above, a corresponds to the spin of the incident particle 
whereas a is a sign index which arises when solving for the 
integral over 0. Inserting this equation into the expression for 
Jy, one arrives at: J y = j^ 2 dQikpsmQ/m^RsHr^Q)] 2 - 

|r°(0)| 2 + e~ 2l * r £ a aarg(a0)}. This expression vanishes 
identically only when the following condition is satisfied: 

/?(6) = r-°(9) (3) 

IfEq. <|3j is not satisfied, i.e. r° ^ r^ °, the transverse spin- 
current remains finite. We underline here that both Eq. (|2|i and 
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FIG. 2: (Color online) Plot of the transverse spin-current J7 V appear- 
ing in a normal metal | normal metal bilayer due to interface spin-orbit 
scattering. We show results for a barrier width of (a) L = 0. 1 nm and 
(b) L = 2 nm, using in both cases a strong barrier potential Zq = 5. 

Eq. ([3]l are completely general expressions that have been de- 
rived without making any assumption about the barrier region 
and the details of the spin-orbit coupling present there. Let 
us now return to our specific model with a spin-orbit coupled 
barrier above contacted to two normal regions and use the po- 
tential V_. We then find the following specific expressions: 
r° = 5 OT (<7P+ -kip- )/{qp+ +klp y ), p± = 1 ±e 2ik < L (k x - 
q)/(k x + q), y = ±1 =t,4 for the reflection coefficients r® = 
{r°,r°}. Above, kl = ^2m(E F + yU soc sin 6 - Uq) - k] is the 
wavevector in the barrier region with k y = k F sin 0. Here, U soc 
and Uq correspond to the spin-orbit coupling and barrier po- 
tential in the central region, respectively. We now prove ex- 
plicitly that the result for r° dictates that the longitudinal spin- 
current vanishes whereas the transverse spin-current is finite. 
By noting that kJ(Q) = k x ^(— 0), we see that the condition 
r° (0) = r^°(— 0) is satisfied. Since Eq. <|2jl is satisfied, it then 
follows that J x = 0. However, Eq. ([3]) is not satisfied, and 
hence we have J y ^ 0. Thus, we conclude that the presence of 
antisymmetric spin-orbit coupling at the interface due to the 
broken inversion symmetry induces a transverse spin-current 
whereas there is no spin flow across the barrier. The trans- 
verse spin-current and its dependence on the non-magnetic 
and spin-orbit coupling potentials is shown in Fig. [2] with 
X F = l/kf. For a normal metal with electron mass with ap- 
proximately its bare value, m ~ niQ, and a Fermi level of order 
e V, Ef ~ 1 e V, we obtain kp — 5 nm~ 1 . We show results for 
a barrier width of (a) L = 0.1 nm and (b) L = 2 nm, quan- 
tifying the potential in the barrier region by the dimension- 
less parameters Zq — 2mUQ/k F and Z soc — 2mU soc /k F . As 
seen, the spin-current increases in magnitude with Z soc and 
oscillates as a function of distance penetrated into the nor- 
mal metal region. In both (a) and (b), we have set the atomic 
barrier potential to be strong, Zq = 5. The main difference be- 
tween these two cases is that the spin-current decays to zero 
in the latter case whereas it saturates at a constant value in 
the former case. In order to understand this physically, we 
note that an increase of the barrier width L inevitably gener- 
ates less transmission through the bilayer such that the prob- 
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FIG. 3: (Color online) Plot of the transverse spin-current J y appear- 
ing in a normal metal | superconductor bilayer due to interface spin- 
orbit scattering. In (a) and (b), we plot the spin-current at Fermi level 
for Z soc = 0.0 to 0.3 in steps of 0. 1 for Zq = 0.5 and Zq = 5, respec- 
tively. In each panel, the curves increase in magnitude with Z soc . In 
the side-panels, a contour-plot of the spin-current in the x-Z soc plane 
is given for the cases of (d) high Zq = 0.5, (e) intermediate Zq = 1, 
and (e) low Zq = 5 interface transparency. 



ability for reflection approaches unity. Examining the analyt- 
ical expression for the transverse spin-current written above, 
it is clear that the constant terms (x-independent) cancel each 
other in this scenario, whereas the oscillatory term remains 
due to the different phases of the spin-t and spin-| reflection 
coefficients. Eq. Q and Q constitute a general criterion for 
the appearance of a pure transverse spin-current in a hybrid 
structure. Although there is no net spin-current flowing across 
the barrier, the interfacial spin-orbit scattering will invariably 
generate a transverse spin-current flowing parallel to the bar- 
rier according to the expression for %. This spin-current will 
oscillate with the distance |jc| away from the barrier, and in 
general die out on a length-scale dictated by spin-relaxation 
processes due to e.g. magnetic impurities and spin-flip scat- 
tering. 

Let us now consider an alternative model for the spin-orbit 
scattering at the interface region to demonstrate that the ap- 
pearance of a transverse spin-current does not pertain to the 
specific model of a spin-orbit coupled barrier used above. 
Moreover, we turn our attention to the scenario where one 
of the contact materials is superconducting in order to see 
how the process of Andreev reflection influences the trans- 
verse spin-current. We employ the potential of the form Bill : 
V = [Vol + Vsoc" • (G x k)/kf]8(r — ro), where ro is the po- 
sition of the interface, and consider for concreteness a bi- 
layer consisting of a normal metal and an s-wave supercon- 
ductor. This enables the electron-hole conversion by means 
of Andreev-reflection. Employing a scattering matrix for- 
mulation, we write down the incoming waves with spin a 
from the normal metal as: y% = (1,0)^*$* + ]$ (l,0)e -i *8* + 
r%(0, lje****. In the superconductor, the transmitted wave- 
function reads: \\% = t%(<&,a)j , w + t%(a,eP)er ]i « x . We 
incorporate any Fermi-vector mismatch as an effective in- 



crease in Vq. The boundary conditions at the interface lo- 
cation x = allow for a determination of the scattering co- 
efficients {rf ,r°,t° ,t%}. These read: = \\t^ and 3 T (\|/g - 
\\f%) = 2m(V - oV S0C sin8)\|$. Upon defining £>° = (Z a ) 2 - 
e 2i P [1 + (Z°) 2 ] where Z° = (Z / cos 9) - oZ S0C tan 9 and Zj = 
mVj/kf, we obtain the following compact expressions for 
the scattering coefficients: rf = Z°(e 2i P - l)(i + Z°)/£>°, 
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We emphasize that these equations are not simply a spin- 
polarized version of the Blonder-Tinkham-Klapwijk result 
ifTZll . The interface spin-orbit coupling considered here gives 
rise to very different scattering coefficients than in the sce- 
nario of e.g. a ferromagnet | superconductor bilayer or even a 
normal | superconductor layer with a magnetic interface where 
majority and minority spin carriers are shifted by an exchange 
field. 

Incorporating the possibility of having Andreev- 
reflection, one finds that an additional term 
(-1) x / " /2 d8£ F sin8/(m)i; a a|r£(8)| 2 is added to %. 
Whereas this term vanishes when Z soc = 0, it remains finite 
in the presence of spin-orbit coupling and moreover enters 
with opposite sign to the probability for normal reflection in 
the expression for jL Thus, whereas it is well-known that 
Andreev-reflection suppresses the longitudinal spin-current J x 
in equal footing as normal reflection without branch-crossing, 
the Andreev-reflection probability enters with opposite sign 
to the normal reflection for the transverse spin-current. The 
physical reason for this can be understood in terms of the 
retroreflective nature of the Andreev scattering process in 
the present system: the back-scattered hole has a group 
velocity opposite to its momentum, and hence propagates in 
the opposite direction of a normally reflected electron along 
the interface. Interestingly, it immediately follows that the 
Andreev-reflection coefficient in undoped graphene, where 
this scattering process can become specular iTHl . would 
contribute to the transverse spin-current in a different way 
than the present case. Thus, the transverse spin-current is ac- 
tually sensitive to the exact nature of the Andreev-scattering 
process, which might be an interesting venue to explore 
further. 

We are now in a position to again prove the existence of 
a non-zero transverse spin-current. The interface spin-orbit 
coupling does not violate time-reversal symmetry, and direct 
inspection of the above result shows that r°(8) = r~ a (—Q) 
holds at the same time as r° (9) 7^ r<r°(9) whenever Z soc 7^ 0. 
Based on our previous analysis, it follows that the transverse 
spin-current J7 V 7^ 0. The inclusion of Andreev reflection al- 
ters the magnitude of the spin-current. We underline here 
that the superscript o for the reflection coefficient denotes 
both the spin of the incident particle and the spin of the re- 
flected particle, such that = r° e , which establishes the cor- 
respondence to Eqs. Q and Q. It is also worth to note 
that the transmission coefficients satsify t® (9) = f~°(— 8) and 



5 (— 9) which can be analytically shown, leading 



to a non-zero spin current on the superconducting side near 
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the interface. 

Our treatment of the generated spin-current so far has been 
quite general as there are only two variational dimensionless 
parameters: the non-magnetic barrier potential Zq and the 
strength of the interface spin-orbit interaction Z soc . We now 
plot the normalized spin-current J7 y to investigate how it de- 
pends on these parameters, and also to see how it decays in- 
side the bulk of the materials as one moves away from the 
interface. In Fig. [3j we consider the behavior of J y in an 
N|S bilayer for the case of a high (a) and low (b) interface 
transparency. In each case, several values of the spin-orbit 
scattering parameter Z soc are considered. As seen, the trans- 
verse spin-current persists in both the N and S region and in- 
creases in magnitude alongside with Z soc . In each of the right 
subpanels of Fig. [3] we provide a contour-plot of the spin- 
current to illustrate its behavior as a function of Z soc and po- 
sition x simultaneously. As seen, the spin-current experiences 
a sharp increase right around Z soc ~ Zq, in effect where the 
spin-orbit scattering potential outgrows the non-magnetic bar- 
rier potential. We also find (not shown) that the spin-current is 
enhanced by the Andreev-reflection process compared to the 
case where the system is non-superconducting, A->0, when 
the interface transparency is high, i.e. low values of Zq. With 
increasing barrier strength Zq, the Andreev-reflection process 
is suppressed and the enhancement vanishes. 

One should bear in mind that detection schemes for spin- 
currents generated by spin-orbit coupling are subject to con- 
straints that arise from time-reversal symmetry in conjunction 
with the Onsager relations fl4l . For instance, it has been 
shown that whereas the detection of spins in a two-probe con- 
ductance setup cannot be accomplished by means of a magne- 
tization switch, the detection of a torque signal generated by 
spin-flow is not prohibited by the symmetry fT51 . 

Finally, we would like to comment on recent works which 
have considered how interfacial spin-orbit scattering arises. 
Whereas structural-inversion asymmetry (SIA) induced spin- 
orbit scattering traditionally has been discussed in the con- 
text of two-dimensional electron gases, ab initio studies have 
also discussed this aspect for metallic surfaces |[T6l and 
ferromagnet| semiconductor junctions (Fe|GeAs) 1171 . More- 
over, a recent experiment demonstrated the spin-orbit scat- 
tering induced spin-torque by SIA in a ferromagnet| insulator 
junction (Co|A10. v |Pt) [18] originally discussed in Q9)- We 
believe that our model for the interfacial scattering should 
have merit since it has been shown that the surface spin-orbit 
scattering on e.g. Au(lll) lfl6l is closely analogous to the 
Rashba interaction in two-dimensional electron gases. Hav- 
ing stated this, we also mention that the effect predicted in 



this Letter does not depend on the explicit Rashba-form used 
for the spin-orbit potential - it only depends on the asymmetry 
present for structural inversion. The analytical calculations 
that we have presented serve as a proof-of-principle for the 
generic effect of a transverse spin-current, whereas a quanti- 
tative estimate of the strength of the interfacial spin-orbit scat- 
tering and the exact plane of the bilayer junction which will 
support the asymmetric crystal-potential would be a task for 
ab initio studies. 

In summary, we have theoretically demonstrated the ap- 
pearance of transverse pure spin-current in any hybrid struc- 
ture, due to broken inversion symmetry near the interface. We 
have provided the condition to realize this effect, and calcu- 
lated the spin-current for two bilayer systems. The predicted 
transverse spin-current may be probed either via conversion 
into an electrical signal by the inverse spin-Hall effect l20ll2TI 
or through optical studies such as non-local Faraday rotation 
ll22ll and polarized light-emission [23 1. 
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